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Abstract. A resistance network is a connected graph (G, c) . The conductance 
function Cxy weights the edges, which arc then interpreted as conductors of 
possibly varying strengths. The Dirichlet energy form £ produces a Hilbert 
space structure T-Le on the space of functions of finite energy. 

The relationship between the natural Dirichlet form £ and the discrete 
Laplace operator A on a finite network is given by £{u, v) = {u, Au)2, where 
the latter is the usual l"^ inner product. We describe a reproducing kernel 
{vx} for £ and use it to extend the discrete Gauss-Green identity to infinite 
networks: 

G bd G 

where the latter sum is understood in a limiting sense, analogous to a Riemann 
sum. This formula yields a boundary sum representation for the harmonic 
functions of finite energy. 

Techniques from stocheistic integration allow one to make the boundary 
bdG precise as a measure space, and give a boundary integral representation 
(in a sense analogous to that of Poisson or Martin boundary theory). This 
is done in terms of a Gel'fand triple S C T-Ls ^ S' and gives a probability 
measure P and an isometric embedding of He into L^(S'',P), and yields a 
concrete representation of the boundary as a set of linear functionals on S. 
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1. Introduction 

There are several notions of "boundary" as "points at infinity" associated to infi- 
nite graphs. Some of these come directly from graph theory, like the notion of graph 
ends [PW90,Die06] or ideal boundary [MY89,Yam86]. Others come from by way of 
the associated reversible Markov process, the random walk associated to the graph, 
like topological notion of Martin boundary [PW87,Saw97] or its measure-theoretic 
refinement, the Poisson boundary [KW07, Kai91]. There are also less well-known 
ideas, like the discrete Royden boundary [KY88] and discrete Kuramochi bound- 
ary [MY97]. Interrelations amongst these concepts are detailed in two excellent 
collections of notes: [WoeOO] and [Soa94]. This material has its roots in minimal 
surface theory, probability theory, ergodic theory, and group theory, and the cen- 
tral ideas are often analogues of a corresponding notion for continuous domains 
(manifolds, Lie groups, etc.). 

This paper gives a brief account of a new type of boundary developed in 
[JP09b, JP09c, JP09a] and [JP08] which we call the resistance boundary; it is de- 
noted bd G. It bears many similarities to the Martin and Poisson boundaries, but 
pertains to a different class of functions: the functions of finite energy. Let G be 
a resistance network (i.e., a connected simple weighted graph) with vertex set G*^ 
and edges determined by a symmetric conductance function c which weights the 
edges: c^y = Cy^ > 0, and c^y > iff there is an edge from x to y, which is denoted 
X ^ y. The energy of a function m : G ^ C is then defined to be 



For the most part, it suffices to work with K- valued functions (see Remark 2.14 
in particular). However, we will need C-valued functions for some applications of 
spectral theory in §5. 

Under suitable hypotheses, if ft, is a bounded harmonic function on X, then 
Poisson boundary theory provides a measure space {dX, fj.) with respect to which 
one has an integral representation of h in terms of a kernel k : X x dX — ^ C: 




(1.1) 




(1.2) 
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where h is the extension of h to dX, in some sense. This paper provides a synopsis 
of how one can obtain a similar representation for the harmonic functions of finite 
energy. However, instead of using ergodic theory or (topological) conipactifications, 
we take an entirely different approach: operator theory and functional analysis. 

After embedding the resistance network into a certain Hilbert space, we con- 
struct a space of distributions (i.e., generalized functions) on that Hilbert space. 
We then show that this space of distributions contains the boundary of the orig- 
inal network, in the sense that it supports integral representations of harmonic 
functions on the network. We work with the energy space, a Hilbert space whose 
inner product is given by the sesquilinear form associated to £ by polarizing (1.1): 

{u,v)£ ^ £{u,v) := ^ c^y (^u{x) ~ u{y)^ (v{x) - v{y)Y (1.3) 

We construct a reproducing kernel for this Hilbert space, and then use it to obtain 
a Gel'fand triple 

SsC-HsC S',. (1.4) 

Here, is a £-dense subspace of which is also equipped with a strictly finer 
"test function topology" (defined in terms of the domain of the Laplacian) , and the 
space is the dual space of Se with respect to this finer topology; the specifics are 
discussed further just below. For now, however, let us eschew technical details and 
just say that S'^ is strictly larger than Tig, and it is in S'g that the boundary bdG 
lies. This framework allows us to invoke Minlos' theorem and Wiener's isomet- 
ric embedding theorem, powerful tools from the theory of stochastic integration. 
Boundary theory usually involves an enlargement of the original space, either by 
topological means (e.g., by compactification or completion, in the case of Martin 
boundary) or by measure-theoretical means (e.g., by taking the measurable hull 
of an equivalence relation, as in Poisson boundary). For the resistance boundary 
bd G, we enlarge Tig (the Hilbert space representation of the resistance network) 
by embedding it into S'^ via the inclusion map. 

Definition 1.1. The Laplacian on a resistance network (G, c) is the linear difference 
operator A which acts on a function i; : G" — > C by 

(Ai;)(x) := ^c,^(i;(x)-t-(y)). (1.5) 

A function w : G*' — > C is harmonic iff Av{x) — for each x E G'~'. 

Note that we adopt the (physicists') sign convention in (1.5) (so that the spec- 
trum is nonncgative) and thus our Laplacian is the negative of the one commonly 
found in the PDF literature; e.g., [KigOl, Str06]. 

The study of resistance boundaries begins with the following well-known iden- 
tity for finite networks. 
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Proposition 1.2. Let G he a finite network. For functions u,v on , 

£{u,v) ^ ^ u{x)Av{x). (1.6) 

The right-hand side of (1.6) is often denoted by {u, Av)2- Theorem 3.3 gives 
a broad extension of Proposition 1.2 to a certain domain Ai (see Definition 2.17). 
Extensions of this type have been studied before (see [Mac80, KY89]), but only 
with regard to determining conditions that ensure £{u, v) — {u, At;)2. By contrast, 
we are more interested in the situation for which it is replaced by 

{u,v)£=J2^Av+Y,u§^- (1-7) 

G° bdG 

Theorem 3.3 gives conditions under which (1.7) holds; the notation bdG and ^ 
are explained precisely in Definition 3.1 and Definition 3.2. In particular, (1.7) 
holds for any u € when v lies in a certain dense subspace of which we 
denote by A4. The space A4 was introduced in [JP09b] for this purpose and also 
to serve as a dense domain for the possibly unbounded Laplace operator, which will 
be useful later for the construction of Se- We call (1.7) the discrete Gauss-Green 
identity by analogy with 

/ WuWvdV^- / uAvdV+ / u§^dS. 
Jn Jn Jan ^ 

The space Hs consists of potentials (functions on the vertices of G, modulo 
constants; see Definition 2.5) and enjoys an orthogonal decomposition into the 
subspace J-in of finitely supported functions and the subspace Harm of harmonic 
functions; this is given precisely in Definitions 2.9-2.11 and Theorem 2.12. It turns 
out that has a reproducing kernel {vj:}x£G°- for any u S Jig, one has 

{vx,u)£ = u{x) - u{o), Vx e G°, 

where o G G'^ is a fixed reference point. Since the reproducing kernel behaves 
well with respect to (orthogonal) projections P, we also have reproducing kernels 
{fx}x£G° foi' s-iid {hx}x£G° for 'Harm, where 

fx '■= PMn''^x^ S-Hd hx ■= Pj^armVx- 

In Theorem 5.1, we apply (1.7) to the reproducing kernels {hx}x£G° for 
Harm, and find that for all h G Harm, 

h{x)^h{o)^Y.^'-t- (1-8) 

bdG 

This direct analogue of (1.2) first appeared in [JP09b, Cor. 3.14]. Formula (1.8) 
gives a boundary sum representation of harmonic functions, but the boundary 
sum in (1.8) is understood only as a limit of sums taken over boundaries of finite 
subnetworks. Comparison of (1.8) and (1.2) makes one optimistic that bdG can 
be realized as a measure space which supports a measure corresponding to 
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thus replacing the sum in (1.8) with a integral. In Corollary 5.19, we extend (1.8) 
to such an integral representation for which (1.8) is analogous to a Riemann sum. 

The primary difference between our boundary theory and that of Poisson 
and Martin is rooted in our focus on Hs: both of these classical theories concern 
harmonic functions with growth/decay restrictions. By contrast, provided they 
neither grow too wildly nor oscillate too wildly, elements of may be unbounded 
and may fail to remain nonnegative. From [ALP99], it is known that functions 
which are £-limits of finitely supported functions must vanish at oo (except for a 
set of measure with respect to the usual path-space measure); however see [JP08, 
Ex. 13.10] for an unbounded harmonic function of finite energy. Note, however, 
that functions of finite energy can always be approximated in by bounded 
functions; cf. [Soa94, §3.7]. 

Just as for Martin and Poisson boundaries, the resistance boundary essen- 
tially consists of different limiting behaviors of the (transient) random walk on 
the network, as the walker tends to infinity. It turns out that recurrent networks 
have no resistance boundary, and transient networks with no nontrivial harmonic 
functions have exactly one boundary point (corresponding to the fact that the 
monopole at x is unique; see Definition 2.15). In particular, the integer lattices 
(Z*^, 1) each have 1 boundary point for d > 3 and boundary points for d = 1,2. 
Further examples are discussed in §6. 

Outline. §2 recalls basic definitions and some previously obtained results. In par- 
ticular, we give precise definitions for the Laplace operator A, the energy space "H^, 
the reproducing kernel {vx}, monopoles Wx, the monopolar domain Ai, and we 
discuss the Royden decomposition of Us into the finitely supported functions and 
the harmonic functions. §3 states the discrete Gauss-Green identity and gives the 
definition of the boundary sum X^bdG ^ ^ limit of sums. Some implications 
of the discrete Gauss-Green identity are given, including several characterizations 
of transience of the random walk on the network. §4 gives the definition of effec- 
tive resistance, and discusses how this metric can be extended to infinite networks 
in different ways the free resistance R^{x,y) and wired resistance R^{x,y). §5 
discusses the boundary sum representation for elements of Harm as introduced in 
(1.8). This section also gives an overview of the theory of Gel'fand triples, Minlos' 
theorem, and Wiener's theorem, and how these enable one to obtain a Gaussian 
probability measure on the space S'^ alluded to in (1.4). §5 gives the boundary 
integral representation of elements of Harm: an integral version of (1.8) which is 
an 'Hf-analogue of (1.2). §6 contains several examples which illustrate our results. 

Boundary theory is a well-established subject; the deep connections between 
harmonic analysis, probability, and potential theory have led to several notions of 
boundary and we will not attempt to give complete references. However, we rec- 
ommend [Saw97] for introductory material on Martin boundary and [WoeOO] for 
a more detailed discussion. Introductory material on resistance networks may be 
found in [DS84] and [LPIO], and [Kig03] gives a detailed investigation of resistance 
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forms (a potential-theoretic generalization of resistance networks). More specific 
background appears in [Lyo83, Car73] and the foundational paper [NW59]. With 
regard to infinite graphs and finite-energy functions, see [Soa94, WoeOO, SW91, 
CW92,Dod06,PW90,PW88,Woc86,Tho90]. Applications to analysis on fractals 
can be found in [KigOl, Str06]. For papers studying fractals as boundaries of net- 
works or Markov processes, see [DS99,DS01,DS02,Kai03,LW09, JLW10,Kig09]. 

2. The energy space 1-Ls 

We now proceed to introduce the key notions used throughout this paper: resis- 
tance networks, the energy form E, the Laplace operator A, the energy space 
the reproducing kernel {vx}, and their elementary properties. 

Deflnition 2.1. A resistance network is a connected graph (G, c), where G is a graph 
with vertex set , and c is the conductance function which defines adjacency by 
a; ^ y iff Cxy > 0, for x,y G G^. We assume Cxy = Cyx G [0, oo), and write 
c{x) :— J2yr^x ^xy Wc require c{x) < oo but c{x) need not be a bounded function 
on G", and note that vertices of infinite degree are allowed. The notation c may be 
used to indicate the multiplication operator (cv){x) := c{x)v{x), i.e., the diagonal 
matrix with entries c{x) with respect to the (vector space) basis {Sx}- 

As the letters x, y, z always denote vertices, it causes no confusion to write 
x,y,z£G instead of x,y, z G G". Similarly, u and v will always denote functions 
which map vertices to scalars, e.g., u : G*' — ?> C. 

In Definition 2.1, "connected" means simply that for any x,y £ G, there is 
a finite sequence {xi}2^Q with x = Xq, y = Xn, and Cxi_ixi > 0, i = l,...,n. 
Conductance is the reciprocal of resistance, so one can think of (G, c) as a network 
of nodes G" connected by resistors of resistance c~J^. We may assume there is at 
most one edge from x to y, as two conductors c^^ and c^y connected in parallel 
can be replaced by a single conductor with conductance Cxy = c^y + c^y Also, we 
assume Cxx = so that no vertex has a loop, as electric current will never flow 
along a conductor connecting a node to itself. 

Definition 2.2. An exhaustion of G is an increasing sequence of finite and connected 
subgraphs {Gfc}^^, so that Gk C Gk+i and G = IJ Gfc. Since any vertex or edge 
is eventually contained in some Gk, there is no loss of generality in assuming they 
are contained in Gi, for the purposes of a specific computation. 

Deflnition 2.3. The notation 



is used whenever the limit is independent of the choice of exhaustion {Gk} of G. 
This is clearly justified, for example, whenever the sum has only finitely many 
nonzero terms, or is absolutely convergent as in the definition of £ just below. 




(2.1) 



x&Gk 
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Deflnition 2.4. The energy of functions m, w : — C is given by the (closed, 
bihnear) Dirichlet form 

£{u,v) -^^^Yl - '^iv))i^i^) - '"(y))^ (2-2) 

xeG yeG 

with the energy of u given by £{u) := £{u,u). The domain of the energy is 

domf = {w : G° ^ C ; £(u) < oo}. (2.3) 

Since Cxy = Cyx and Cxy = for nonadjacent vertices, the initial factor of ^ 
in (2.2) implies there is exactly one term in the sum for each edge in the network. 

Deflnition 2.5. Let 1 denote the constant function with value 1 and recall that 
ker£ = CI. The energy form £ is symmetric and positive definite on dom£. Then 
dom£/Cl is a vector space with inner product and corresponding norm given by 

{u,v)£ ■.^£{u,v) and \\u\\£ £{u,uy/^ . (2.4) 

The energy Hilbert space is defined to be 

dom£ dom£ 
- ke^ ^ -CT- (2-^) 

Thus, He consists of potentials: we are not interested in values u{x) as much 
as differences u(x) — u{y). In other words, if u and v are both elements of dom£ 
and there is some constant fc G C such that u{x) ~ v{x) — k for all x E G, then u 
and V are both representatives of the same element (equivalence class) of . 

Deflnition 2.6. Let be defined to be the unique element of Hs for which 

{vx,u)£ — u{x) — u{o), for every u e %£. (2-6) 

The collection {vx}xeG forms a reproducing kernel for T-Lg] cf. [JP09b, Cor. 2.7]. 
We call it the energy kernel and (2.6) shows its span is dense in "Hf. Note that Vo 
corresponds to a constant function, since (uq, u)^ = for every u G %£. Therefore, 
Vo (or o) may often ignored or omitted. 

Definition 2.7. A dipole is any v G T-Lg satisfying the pointwise identity Aw = 5x—5y 
for some vertices x,y E G. The elements of the energy kernel are all dipoles: one 
can check that Avx — 6x ^ 6o as in [JP09b, Lemma 2.13]. 

Remark 2.8. To minimize cumbersome notation, let {x G G} be the default index 
set from now on. That is, we use {fa:} to denote the energy kernel {vx}xeGi ^-nd 
si>an{vx} to denote the set of all linear combinations of elements of {vx}, etc. 

2.1. The finitely-supported functions and the harmonic functions 

Deflnition 2.9. For v G He, one says that v has finite support iff there is a finite 
set F C for which v{x) = fc G C for all a: ^ F, i.e., the set of functions of finite 
support in Jig is spanjJa,}, where Sx is the Dirac mass at x, i.e., the element of 
Us containing the characteristic function of the singleton {x}. Define J-in to be 
the closure of spanj^a;} with respect to £. 
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Remark 2.10. The usual candidate for an orthonormal basis (onb) in ^^(G") would 
be the collection of Dirac masses {Sx}', however, this is not an onb in Hg. One can 
compute from (2.2) that 

{Sx,Sy)£ = £{Sx,Sy) ^ -Cxy, for a; 7^ y, and £{Sx) = c{x), (2.7) 

so that Sx X- with respect to £. Moreover, Theorem 2.12 shows that {Sx} is, in 
general, not even dense in %£. It is immediate from (2.7) that 6x G T-Ls- 

Definition 2.11. The harmonic subspace of T-Lg is denoted 

Uarm {v G Us \ Aw(a;) = 0, for ah x e G}. (2.8) 

Note that this is independent of choice of representative for v in virtue of (1.5). 

The following result is sometimes called the "Royden Decomposition" since 
[Yam79, Thm. 4.1], in reference to Royden's analogous result for Riemann surfaces; 
see [Soa94, §VI], [LPIO, §9.3]^ It follows immediately from [JP09b, Lemma 2.11], 
which states that {Sx,u)£ = Au{x) for any x G G; cf. [JP09b, Thm. 2.15]. 

Theorem 2.12 (Royden decomposition). Hs = J-in ® Harm. 

Definition 2.13. Let fx = PnnVx denote the image of Vx under the (orthogonal) 
projection to J-in. Similarly, let hx = Pnarm'^x denote the image of Vx under the 
projection to Harm. 

Remark 2.14 (Reproducing kernels for J-in and Harm). The reproducing kernel 
property behaves well with respect to orthogonal projections, and consequently, 
{fx} is a reproducing kernel for J^in, and {hx} is a reproducing kernel for Harm. 
While we will need complex-valued functions for some results obtained via spec- 
tral theory, it will usually suffice to consider M-valued functions because the re- 
producing kernels elements Vx,fx,hx all have M- valued representatives [JP09b, 
Lemma 2.24]. 

2.2. Monopoles 

Definition 2.15. A monopole at x E G is an element Wx G Hg which satisfies 
Awx{y) — Sxy, where Sxy is Kronecker's delta. In case the network supports 
monopoles (that is, if the above Dirichlet equation admits finite-energy solutions), 
let Wo always denote the unique energy-minimizing monopole at the origin. 

With Vx and fx = PnnVx as above, we indicate the distinguished monopoles 

wl := Vx -\- Wo and w^x ■= fx + Wo- (2.9) 

Remark 2.16. Note that Wo G J-in, whenever it is present in Hs, and similarly 
that wlc is the energy-minimizing monopole at x. To see this, suppose Wx is any 
monopole at x. Since Wx G Hg, write Wx — f + h hy Theorem 2.12, and get 
£iwx) = £{f) -\-£{h). Projecting away the harmonic component will not affect the 
monopole property, so w4 = Pj=inWx is the unique monopole of minimal energy. 



'^This name is also sometimes associated with the corresponding (nonorthogonal) decomposition 
for the "grounded energy form"; see Remark 3.11. 
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The Green function is g{x, y) — Wy{x), where w" is the representative of Wy which 
vanishes at oo. 

Definition 2.17. The dense subspace oiHs spanned by monopoles (and dipoles) is 
M spaii{vx} + span{w^, (2-10) 
Let Am be the closure of the Laplacian when taken to have the dense domain Ai . 

Since A agrees with Am pointwise, we may suppress reference to the domain 
for ease of notation. It is shown in [JP09b, Lemma 3.5] that Avi is Hermitian with 
{u, Aj\^u)£ > for ah u G M.. When given a pointwise identity Au = w, there is 
an associated identity in He, but one must use the adjoint: Au{x) — v{x) for ah 
a; G G if and only \i v = A^u in Hs [JP09b, Lemma 3.7]. Note that Am may 
have defect vectors; such an object is an element of the Hilbert space Jig (though 
clearly not an element of A^) which has a representative u satisfying 

Amu{x) = —u{x),\fx e G, and u € domA^. 

See [JP09d, §4.2] or [JP08, §13.4]. While it is always the case that the (possibly 
unbounded) operator Am is Hermitian (i.e. Am ^ Am)' ^^^^ shows that Am may 
fail to be self-adjoint (i.e. Am ~ A^). 

Remark 2.18 (Monopoles and transience). The presence of monopoles in Hs is 
equivalent to the transience of the simple random walk on the network with tran- 
sition probabilities p{x,y) = Cxy/c{x): note that if Wx is a monopole, then the 
current induced by Wx is a unit flow to infinity with finite energy. It was proved 
in [Lyo83] that the network is transient if and only if there exists a unit current flow 
to infinity; see also [LPIO, Thm. 2.10]. Moreover, it is shown in [JP09b, Lemma 3.6] 
that when the network is transient, Ai contains the spaces sp&n{vx}, span{/a;}, and 
span{/ia;}, where fx = Pjin^x and hx = PnarmVx- When Harm — (in particular, 
when the network is not transient), fx — Vx and so = spanju^:} = span{/a;} 
trivially. 



3. The discrete Gauss-Green formula 

In Theorem 3.3, we establish a discrete version of the Gauss-Green formula which 
extends Proposition 1.2 to the case of infinite graphs; the scope of validity of this 
formula is given in terms of the space Ai of Definition 2.15. The appearance of a 
somewhat mysterious boundary term alluded to in (1.7) prompts several questions 
which are discussed in Remark 3.6. 

3.1. Relating A to f 

Definition 3.1. If is a subgraph of G, then the boundary of H is 

hdH -.^ {x e H [3y e H\y ^ x}. (3.1) 
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The interior of a subgraph H consists of the vertices in H whose neighbours also 
he in H: 

iiAH := {x & H \ y ^ X =^ y & H} = H\hdH. (3.2) 
For vertices in the boundary of a subgraph, the normal derivative of v is 

^{x) := J2 c-vi^i^) - ^(y))^ for xehdH. (3.3) 

yeH 

Thus, the normal derivative of v is computed like Av{x), except that the sum 
extends only over the neighbours of x which lie in H. 

Definition 3.1 will be used primarily for subgraphs that form an exhaustion 
of G, in the sense of Definition 2.2. 

Definition 3.2. A boundary sum is computed in terms of an exhaustion {Gk} by 

bdG bdGfc 

whenever the limit is independent of the choice of exhaustion, as in Definition 2.3. 
Theorem 3.3 (Discrete Gauss-Green Formula). If u ^ Hg and v E M, then 

{u,v)s=Y.uAv+J2ul^- (3-5) 

G bdG 

Corollary 3.4. For all u e domAvi, '^qAu ~ ^SbdG f^' '^^^-"'^ discrete 
Gauss- Green formula (3.5) is independent of choice of representatives. 

Remark 3.5. The proof of Theorem 3.3 follows from taking limits of 

^ uix)Avix)+ ^ uix)§^{x). 

xeGk xGbdGfe 

Thus, the decomposition (3.5) is true for all u,v G Hg, but is meaningless if it 
takes the form oo — oo. A key point of Theorem 3.3 is that for u, v in the specified 
domains, the two sums are both finite and independent of choice of exhaustion. 
However, the specific value of each sum is dependent on the choice of representative 
for u; this motivates Definition 3.9. 

It is also clear that (3.5) remains true much more generally than under the 
specified conditions; certainly the formula holds whenever J^xgG \u{x)Av{x)\ < oo. 
Unfortunately, given any hypotheses more specific than this, the limitless variety 
of infinite networks almost always allows one to construct a counterexample; i.e. 
one cannot give a condition for which the formula is true for all u € "Hf , for all net- 
works. Nonetheless, the formula remains true and even useful in many specific and 
general contexts. For example, it is clearly valid whenever w is a dipole, including 
all those in the energy kernel. We will also see that it holds for the projections of 
Vx to J-in and to 'Harm. Consequently, for v which are limits of elements in AA, 
we can use this result in combination with ad hoc arguments. 

A formula similar to (3.5) appears in [DK88, Prop 1.3]; however, these au- 
thors apparently do not pursue the extension of this formula to infinite networks. 
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Another similar result appears in [KY89, Thm. 4.1], where the authors give some 
conditions under which (1.6) extends to infinite networks. The main differences 
here are that the scope of Kayano and Yamasaki's theorem is limited to a subset 
of what we call J-in^ and that Kayano and Yamasaki are interested in when the 
boundary term vanishes; we are more interested in when it is finite and nonvanish- 
ing; see Theorem 3.10, for example. Since Kayano and Yamasaki do not discuss the 
structure of the space of functions they consider, it is not clear how large the scope 
of their result is; their result requires the hypothesis X^xsg |w(x)Au(x)| < oo, but 
it is not so clear what functions satisfy this. By contrast, we develop a dense sub- 
space of functions on which to apply the formula. Furthermore, in the forthcoming 
paper [JP09c], we show that these functions arc relatively easy to compute. 

Remark 3.6. We refer to X^bdG^fn "boundary term" by analogy with 

classical PDE theory. This terminology should not be confused with the notion of 
boundary that arises in the discussion of the discrete Dirichlet problem, where the 
boundary is a prescribed subset of . As the boundary term may be difficult to 
contend with, it is extremely useful to know when it vanishes, for example: 

(i) when the network is recurrent (Theorem 3.10), 

(ii) when v is an element of the energy kernel [JP09b, Lemma 5.8], 

(iii) when M,ti, Au, Aw lie in £^ [JP09b, Lemma 5.12], and 

(iv) when either u or d has finite support [JP09b, Lemma 5.16]. 

3.2. More about monopoles and the space M 

This section studies the role of the monopoles with regard to the boundary term of 
Theorem 3.3, and provides several characterizations of transience of the network, 
in terms the operator-theoretic properties of A_a4 . 

Note that if /i e Harm satisfies the hypotheses of Theorem 3.3, then £{h) = 
SbdG^is- other hand, £{u) = uAit for all u G Hs iff the network 

is recurrent, as stated in Theorem 3.10. With respect to Hs = Fin ® Harm, 
this shows that the energy of finitely supported functions comes from the sum 
over G, and the energy of harmonic functions comes from the boundary sum. 
However, for a monopole Wx, the representative specified by Wx{x) = satisfies 
^('"^) = Sbdc'^^is^ representative specified by Wx{x) = £{wx) satisfies 

^('"^) — TliG w^w. Roughly speaking, a monopole is therefore "half of a harmonic 
function" or halfway to being a harmonic function. A further justification for this 
comment is given by Corollary 3.8 (the proof shows that a harmonic function can 
be constructed from two monopoles at the same vertex, see [JP09b, Cor. 4.4]). 
The general theme of this section is the ability of monopoles to "bridge" the finite 
and the harmonic. 

Theorem 3.7 ( [Soa94, Thm. 1.33]). Let u he a nonnegative function on a recurrent 
network. Then u is superharmonic if and only if u is constant. 

It follows from Theorem 3.7 that Harm ^ implies the existence of a mono- 
pole in He, i.e., the transience of the network; cf. [JP09b, Cor. 4.3]. However, it 
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turns out that a nontrivial harmonic function can only exist when there is more 
than one monopole. 

Corollary 3.8. Harm ^ Q iff there are at least two linearly independent monopoles 
at one (equivalently, every) vertex x. 

Definition 3.9. The phrase "t/ie boundary term is nonvanishing" indicates that 
(3.5) holds with nonzero boundary sum when applied to {u,v)£, for every repre- 
sentative of u except one; namely, the one specified by u{x) — {u, 'w^)£. 

Recall from Remark 2.18 that the network is transient iff there are monopoles 
in Hg. From the Discrete Gauss-Green theorem, we obtain three more criteria for 
transience of the random walk. 

Theorem 3.10. The random walk on the network (G, c) with transition probabilities 
p{x, y) = is transient if and only if any of the following equivalent conditions 
are satisfied: 
(i) the boundary term is nonvanishing, 

(a) fk '■= (cfe + ^)~'^5x is weak-* convergent for some sequence £k — >■ 0, or 
(Hi) (ran A;|;,)'='' = Mn. 

Note that on any network, (ranAA^)*^'" C Tin and hence Harm C kerAJ^; 
cf. [JP09b, Lemma 4.8]. 

Remark 3.11. An alternative approach to studying the space of finite-energy func- 
tions comes by considering the grounded inner product 

{u, v)o := u{o)v{o) + (u, v)£, 

which makes dom £ into a Hilbert space D which we call the grounded energy space. 
This approach is discussed in [LPIO], [Soa94] and [KY89,KY84,MYY94,Yam79]. 

Let Dp be the closure of spanj^j,} in D. If is the projection to Dq, and 
it is applied to the constant function 1 then 1 = 1 if and only if the network is 
recurrent. In fact, when the network is transient, then (modulo additive constants) 
both Pogl and Pd„1 are scalar multiples of monopoles at o. The space Dq is 
spanned by monopoles and harmonic functions. See [JP09b, §4.1] for more details. 

4. Effective resistance 

There is a natural notion of distance on finite networks, which is defined in terms 
of resistance. Consider each edge of the network to be an electrical resistor of 
resistance c~y. The effective resistance metric R{x,y) is the voltage drop between 
the vertices x and y if a current of one amp is inserted into the network at x 
and withdrawn at y. It is a bit surprising that this actually gives a metric, and 
there are several other equivalent formulations, most of which are well-known. The 
essential reference for effective resistance is [Kig03], but the reader may also find 
the excellent treatments in [Soa94] and [LPIO] to be helpful. 
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Figure 1. Effective resistance as network reduction to a trivial network. 
This basic example uses parallel reduction followed by series reduction; see 
Remark 4.2. 

Theorem 4.1. The resistance R{x, y) has the following equivalent formulations: 

R{x, y) = {v{x) - v{v) : ^v = 5.,~ 5,} (4.1) 

= {£{v) \Av^S,~ Sy} (4.2) 

= 1/ niin{£(w) ; v{x) — 1, v{y) — 0,v £ domS} (4.3) 

= min{K > : \v{x) - v{y)f < k£{v), v G dom8} (4.4) 

= sup{|u(a;) - v{y)\'^ ; £{v) < 1, u G dom£}. (4.5) 

Remark 4.2 (Resistance distance via network reduction). Let G be a finite planar 
network and pick any x,y G G^. Then G may be reduced to a trivial network 
consisting only of these two vertices and a single edge between them via the use 
of three basic transformations: (i) series reduction, (ii) parallel reduction, and (iii) 
the V-Y transform [Epi66, Tru89]. The effective resistance between x and y may 
be interpreted as the resistance of the resulting single edge; see Figure 1. See 
also [KigOl] or [Str06] for the V-Y transform. 

4.1. Resistance metric on infinite networks 

There are challenges in extending the notion of effective resistance to infinite net- 
works. The existence of nonconstant harmonic functions h € domS implies the 
nonuniqueness of solutions to Am = / in Hs, and hence (4.1) and (4.2) are no 
longer well-defined. This issue is studied in detail in [JP09c], and in [Kig03] (by 
very different methods). There are also accounts in [LPIO] and the literature on 
"uniqueness of currents" in infinite networks, e.g. [SW91,Tho90]. 

Two natural choices for extension lead to the free resistance R^ and the 
wired resistance R^ . In general, one has R^{x,y) > R^{x,y) with equality iff 
Harm = 0. Both of these correspond to the selection of certain solutions to Am = 
Sx — (in fact, these can be interpreted as Neumann and Dirichlet boundary 
conditions, respectively; see [JP09c, Rem. 2.23]). Also, both are given in terms of 
limits computed with respect to certain networks associated to an exhaustion, in 
the sense of Definition 2.2. The notation {Gk}'^^ always denotes an exhaustion 
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of the infinite network (G, c), as in Definition 2.2. Since x and y are contained in 
all but finitely many Gk, we may always assume that x,y Cz Gk, V/c. 

Definition 4.3. If is a finite subnetwork of G which contains x and y, define 

Rh{x, y) to be the resistance distance from x to y as computed within H. In other 
words, compute RH{x,y) by any of the equivalent formulas of Theorem 4.1, but 
extremizing over only those functions whose support is contained in H. 

Definition 4.4. Let iJ" C G". Then the full subnetwork on has all the edges 
of G for which both cndpoints lie in H^, with the same conductances. That is, 

4.1.1. Free resistance. 

Definition 4.5. For any subset H'-' C G°, the free subnetwork is just the full 
subnetwork with vertices . That is, all edges of G with endpoints in are 
edges of if^, with the same conductances. Thus, we will denote by H to 
reduce notation. Let Rnix^y) denote the effective resistance between x and y as 
computed in H, as in Definition 4.3. The free resistance between x and y is defined 
to be 

R^ix,y) := lim RG,{x,y), (4.6) 

k^oo 

where {Gk} is any exhaustion of G. 

The name "free" comes from the fact that this formulation is free of any 
boundary conditions or considerations of the complements of the Gk; see [LPIO, 
§9]. Theorem 4.6 is the free extension of Theorem 4.1 to infinite networks. 

Theorem 4.6 ( [JP09c, Thm. 2.14]). For an infinite network G, the free resistance 
R^{x,y) has the following equivalent formulations: 

R^{x,y) = v{x) - v{y), v^v^-Vy (4.7) 

= £{v), v = Vx-Vy (4.8) 

= min{i:>(/) ;• / e Tix, y) andl ^7^,1 (4-9) 

= {mm{£{u) : u e Hf, \u{x) - u{y)\ = 1})"^ (4.10) 

= inf{K > ;■ \v{x) ~ w(y)p < K£{v)yv e "Hf} (4.11) 

= sup{\vix) - v[y)\^ we-Hs, Ms < 1} (4.12) 

Fix x,y G G and define the operator L^y on Hs by L^yV := v{x) ~ v{y). 
Then (4.11)-(4.12) are equivalent to R^{x,y) = \\Lxy\\. 

4.1.2. Wired resistance. 

Definition 4.7. Given a finite full subnetwork H of G, define the wired subnetwork 
by identifying all vertices in G" \ to a single, new vertex labeled oo. Thus, 
the vertex set of is iJ" U {oo^f}, and the edge set of includes all the 
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Figure 2. Comparison of free and wired exhaustions for the example of 
the binary tree; see Definition 4.5 and Definition 4.7. Here, the vertices of Gk 
are all those which lie within k edges ( "steps" ) of the origin. If the edges of G 
all have conductance 1, then so do all the edges of each G^ and G^ , except 
for the edges incident upon oo^ = coq^. , which have conductance 2. 



edges of H, with the same conductances. However, if a; e i/" has a neighbour 
y G \ then also includes an edge from x to oo with conductance 

The identification of vertices in G\ may result in parallel edges; then (4.13) cor- 
responds to replacing these parallel edges by a single edge according to the usual 
formula for resistors in parallel. 

Let RHw(x,y) denote the effective resistance between x and y as computed 
in iJ^, as in Definition 4.3. The wired resistance is then defined to be 

lim RGf{x,y), (4.14) 

k^OQ ^ 

where {Gfc} is any exhaustion of G. 

The wired subnetwork is equivalently obtained by "shorting together" all 
vertices of H'^, and hence it follows from Rayleigh's monotonicity principle that 
R^{x,y) < R^ix,y); cf. [DS84, §1.4] or [LPIO, §2.4]. 

Theorem 4.8 ( [JP09c, Thm. 2.20]). The wired resistance may be computed by any 
of the following equivalent formulations: 

R'^{x,y)^f{x)-f{y), f = fx - fy (4.15) 

= f(/), f^fx~fy (4.16) 
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= {min{£{v) : \v{x) - v{y)\ = 1, w G Jin})"' (4.17) 

= inf{K > ; \v{x) - u(y)p < k£(u),Vu £ Tin} (4.18) 

sup{|u(a;) - u(y)p : v e Jin, < 1} (4.19) 
Note that (4.15) and (4.16) are equivalent to 

{x, y) — min{v{x) — v{y) \ Au — 6^ ~ 5y^v <E dom£} (4-15') 

= mva{E{v) ; ls.v ^ 5x ~ Sy,v G dom£}. (4.16') 
4.2. von Neumann construction of the energy space Hs 



Let R = 01 R = i?^ . The discussion of the effective resistance is important in 
this paper in two respects. 

(i) Theorem 4.11 shows that Hs is the natural Hilbert space for studying the 
metric space {G,R). 

(ii) The function R{x, y) allows us to construct a probability measure in Theo- 
rem 5.15. 

Both of these results stem from the fact that (free or wired) effective resistance is 
a negative semidefinite function on x G^, as is shown in [JP09c, Thm. 5.4]. 

Definition 4.9. A function AI : X x X ^ M. is negative semidefinite iff for any 
/ : X — >• K satisfying J2xex /(^) ^ ™^ 

f{x)M{x,y)f{y)<0, (4.20) 

x,y^F 

where F is any finite subset of X. 

One can think of M as a matrix and (4.20) as matrix multiplication, von 
Neumann and Schoenberg [vN32,Ber96,Sch38b,Sch38a,BCR84] showed that (4.20) 
is precisely the condition that allows one to embed a metric space into a Hilbert 
space. This theorem also has a form of uniqueness which may be thought of as a 
universal property. 

Theorem 4.10 (von Neumann). Suppose {X, d) is a metric space. There exists a 
Hilbert space % and an embedding w : {X, d) — >■ "H sending x i-^ Wx and satisfying 

d{x,y) ^ \\wx - WyW-H (4.21) 

if and only if d^ is negative semidefinite. 

Furthermore, if there is another Hilbert space K, and an embedding k :% ^ K,, 
with \\kx — kyWx, = d{x,y) and {kx}xex dense in K,, then there exists a unique 
unitary isomorphism U : "H — > /C. 

Theorem 4.11 ( [JP09c, Thm. 5.4]). {G,R^) may be isometrically embedded in 
a Hilbert space, and this Hilbert space is unitarily equivalent to %£. Under this 
embedding, x is mapped to the energy kernel element Vx ■ 

Moreover, {G, i?^) may be isometrically embedded in a Hilbert space, and this 
Hilbert space is unitarily equivalent to J-in. Under this embedding, x is mapped to 

fx ^Tin^x- 
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By this theorem, we see that Hg is the natural choice of Hilbert space for 
studying the metric spaces {G,R^) and {G,R^). 

5. A boundary integral representation for the harmonic functions 

We are motivated by the following result, which follows readily from Theorem 3.3 
and may be found in [JP09b, Cor. 3.14]. 

Theorem 5.1 (Boundary representation of harmonic functions). Foru G spa.n{hx}, 

^a;) = ^«f^+«(o). (5.1) 

bdG 

Proof. Note that u{x) - u{o) = {vx,u)£ = {u,v^)£ = Y.hda'^^ by (2.6). □ 

Formula (5.1) begs comparison with the Poisson integral formula. Recall the 
classical result of Poisson that gives a kernel fc : x dVl M from which a bounded 
harmonic function can be given via 

u{x)^ [ u{y)k{x,dy), y e dfl. (5.2) 
Jan 

One would like to obtain a (probability) measure space to serve as the bound- 
ary of G. We use some techniques from the theory of stochastic integration for 
which it was shown in [Ncl64] that a Hilbert space does not suffice; see [HidSO, §3.1]. 
The workaround is to build a Gel'fand triple 5' C "H C S" (a more precise defi- 
nition appears just below), and construct a suitable probability measure on S". 
In §5.1, we briefly describe the general theory of Gel'fand triples as they apply 
in the current context. In §5.2, we use A to construct a Gel'fand triple for Hg. 
Then in §5.3, we apply the general theory to the Gel'fand triple Se C He C S'g 
and obtain a Gaussian probability measure P on 5g , and an isometric embedding 
Hs ^ L^(5£,P). This allows us to study the boundary bdG as a subset of S'^. 
For the general theory of analysis in Hilbert space, see [Gro67, Gro70] . 

5.1. Gel'fand triples and duality 

In a little more detail, a Gel'fand triple (also called a rigged Hilbert space) is 

ScncS', (5.3) 

where S is dense in H and S' is the dual of S. While 5 is a dense subspace of H 
with respect to the Hilbert norm, it also comes equipped with a strictly finer "test 
function" topology, and it is required that the inclusion mapping of 5* into H is 
continuous with respect to these topologies. Therefore, when the dual S' is taken 
with respect to this finer topology, one obtains a strict containment Jis C 5". It 
turns out that S" is large enough to support a (Gaussian!) probability measure. 

We will give a "test function topology" as a Frechet topology defined via a 
specific sequence of seminorms. It was Gel'fand's idea to formalize this construc- 
tion abstractly using a system of nuclearity axioms [GMS58,Min58,Min59]. This 
presentation is adapted from quantum mechanics. 
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Remark 5.2 (Tempered distributions and the Laplacian). There is a concrete sit- 
uation when the Gel'fand triple construction is especiaUy natural: H = L'^{M.,dx) 
and S is the Schwartz space of functions of rapid decay. That is, each / g S" is a 
C°° smooth function which decays (along with all its derivatives) faster than any 
polynomial as x ±00. In this case, S' is the space of tempered distributions and 
the seminorms defining the Frechet topology on S are 

p™(/) :=sup{|a;V^"H2^)l x e R,0 < k,n < m}, m = 0,l,2,..., 

where J*^"-' is the 71*^^ derivative of /. Then S" is the dual of S with respect to this 
Frechet topology. One can equivalently express S as 

S:={f e L^{R) ; {P^+Q^y'f e L\R)yn}, (5.4) 

where P : f{x) j-j^ and Q : f{x) 1— > xf{x) are Heisenberg's operators. The 
operator + is often called the quantum mechanical Hamiltonian, but some 
others (e.g., Hida, Gross) would call it a Laplacian, and this perspective tightens 
the analogy with the present context. In this sense, (5.4) could be rewritten S := 
domA°°; compare to (5.10) just below. 

The duality between S and S" allows for the extension of the inner product 
on H to a pairing of S and S": 

{■,-)-H-'Hxn^C to S X S' ^R. (5.5) 

In other words, one obtains a Fourier-type duality restricted to S. 

The proof of Theorem 5.15 will require Minlos' generalization of Bochner's 
theorem from [Min63, Sch73]. This important result states that a cylindrical mea- 
sure on the dual of a nuclear space is a Radon measure iff its Fourier transform 
is continuous. In this context, however, the notion of Fourier transform is infinite- 
dimensional; cf. [Lee96]. 

Theorem 5.3 (Minlos). Given a Gel'fand triple 5* C "H C S", there is a bijec- 
tive correspondence between the positive definite functions f on S and the Radon 
probability measures on S' , determined uniquely by the identity 

f{s) = / e*<^^«>« dFfiO, Vs e S, (5.6) 
Js' 

where (•, ■)-f^ is the extended pairing on S x S' as in (5.5). 

Formula (5.6) may be interpreted as defining the Fourier transform of P. We 
apply Minlos' theorem in the standard manner for white noise constructions, and 
obtain the following corollary. 

Corollary 5.4 (White noise). Given a Gel'fand triple S* C "H C S", there is a 
probability measure P on 5" satisfying 

g-i(s,s)„ ^ f gi(^,?>^ ^(^). (5.7) 
Js' 
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In the proof of Theorem 5.15, we show that P in (5.7) is actually a Gaussian 
measure on S'g. The function on the left-hand side of (5.7) plays a special role 
in stochastic integration, and its use in quantization. To see that it is a positive 
definite function on S, we appeal to a famous result of Schoenberg which may be 
found in [BCR84,SW49]. 

Theorem 5.5 (Schoenberg). Let X be a set and let Q : X x X ^ M. be a function. 
Then the following are equivalent. 

1. Q is negative semidefinite. 

2. yt e M+, the function pt{x,y) := e^*'5(2::j/) jg positive definite on X x X . 

3. There exists a Hilbert space % and a function f : X ^ % such that 

Q{x,y) = \\f{x)-f{y)\\l. 

In the proof of Theorem 5.15, we apply Schoenberg's Theorem with t = ^ to 
the resistance metric in the form 

R^{x,y)^\\v^-Vy\\l, (5.8) 

which appears in [JP09c, Thm. 2.13]. Recall from §4.2 that (5.8) is negative semi- 
definite. 

5.2. A Gel'fand triple for Hg 

To apply Minlos' Theorem, we first need to construct a Gel'fand triple for %£] 
we begin by identifying a certain subspace of = dom Aa/i (as given in Defini- 
tion 2.17) to act as the space of test functions, which we denote Ss- 

Definition 5.6. Let Aa^ be a self-adjoint extension of Aa^; since is Hermitian 
and commutes with conjugation (since c is M- valued), a theorem of von Neumann's 
states that such an extension exists. 

Let A^u :— (Atm A,m . . . Am)m be the p-fold product of A^vi applied to m G %£. 
Define dom(A^) inductively by 

dom(AP,) := {u \ A^'^w G dom(AM)}. (5.9) 

Definition 5.7 (Test functions). The (Schwartz) space of potentials of rapid decay 
is 

5£:=dom(A^), (5.10) 
where dom(A^) dom(A^) consists of all u £ Tig for which A^u € %£ 

for any p. 

Definition 5.8 (Distributions). For each p e N, there is a seminorm on Ss defined 

by 

||u||p:=||A^«||£. (5.11) 
Since (dom A^, \\ ■ \\p) is a Hilbert space for each p S N, the system of seminorms 
^ — {II ■ llplpeN defines a Frechet topology on The space 5^ of Schwartz 
distributions or tempered distributions is the (dual) space of P-continuous linear 
functionals on iSg. 
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Remark 5.9. If deg(a;) is finite for each x G G*^, or if ||c|| < oo, then one has Vx S Sg. 
In the first case, this can be proved from the identity Sx = c{x)vx — <^xyVy 
which is given in [JP09b, Lem. 2.22]. In the second case, the bound on c implies 
Am is bounded and hence everywhere-defined. 

When Se contains {vx}, it should be noted that sj>an{vx} is dense in Se 
with respect to £, but not with respect to the Frechet topology induced by the 
seminorms (5.11), nor with respect to the graph norm. One has the inclusions 



Vx 
AmVx 



C 



s 



C 



u 



(5.12) 



where s € 5^ and u e Hg. The second inclusion is dense but the first is not. 
Remark 5.10. Note that Sf and Si- consist of M-valued functions. This technical 



detail is important because we do not expect the integral Jg, e' 



' from (5.6) 



to converge unless it is certain that {u, •) is M-valued. This is the reason for the 
last conclusion of Theorem 5.13. 

Deflnition 5.11. Let X[a denote the usual indicator function of the interval [a, b] C 
M, and let © be the spectral transform in the spectral representation of Am, and let 
E be the associated projection- valued measure. Then define E^ to be the spectral 
truncation operator acting on by 



E„u := &*X,i „i6u = / E{dt) 

Jl/n 



/r. 

Lemma 5.12. With respect to £, Sg is a dense analytic subspace ofHg. 

Proof. This essentially follows immediately once it is clear that £"„ maps Hg into 
iSg. For u £ He, and for any p ^ 1,2,..., 



X''P\\E{dX)u\\l < n^P\\u\\ 



l/r. 



(5.13) 

□ 



So EnU e iSf. It follows that ||u — Enu\\£ — > by standard spectral theory. 

Theorem 5.13. Sg C "He C S'^ is a Geljand triple, and the energy form 
extends to a pairing on x defined by 

(«,Ow:=(A>,A-fO£, (5.14) 

where p is any integer such that \^{u)\ < K\\APu\\£ for all u Cz Sg. This pairing 
on Ss X iSg is equivalently given by 

{■u,Ow= lim f(£;„u), (5.15) 

n— f oo 

where the limit is taken in the topology of S'^ . 

Corollary 5.14. En extends to a mapping En : S'^ — Us defined via {u,En^)£ '■= 
^{E„u). Thus, we have a pointwise extension of (• , •)w to %£ x 5^ given by 

{u,Ow^ \im {u,En£.)£- (5.16) 
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5.3. The Wiener embedding and the space 5^ 

With Theorem 5.13, we have a Gel'fand triple and we are now ready to apply the 
white noise construction of Cor. 5.4. Note that in Theorem 5.15, expectations are 
taken with respect to the variable S S'^, that is, E(/) := J^, /(^)dP(^). 

Theorem 5.15 (Wiener embedding). The Wiener transform W : Hg L'^{S'g,P) 
is given by 

W:v^v, viO:={v,Ow, (5.17) 

and is an isometry. The extended reproducing kernel {wxlajgcf ^ system of 
Gaussian random variables which gives the resistance distance by 

R^{x,y)^E{{i,~Vy)^). (5.18) 

Moreover, for any u, u € "Hf, the energy inner product extends directly as 




(5.19) 



Proof. Since R^{x,y) is negative semidefinite (see [JP09c, Thm. 5.4]), we may 
apply Schoenberg's theorem and deduce that exp(— u|||) is a positive definite 
function on Jig x Jig. Consequently, an application of the Minlos correspondence 
(Theorem 5.3) to the Gel'fand triple established in Lemma 5.12 yields a Gaussian 
probability measure ¥ on S'g. 
Moreover, (5.6) gives 

E(ei(".«>w) = e-5ll«ll£, (5.20) 

provided that ^ is M- valued (so that the integral converges). Therefore, we give 
the proof for the M- valued subspace of iSg (and of iS^), and then complexify in the 
last step via the standard decomposition into real and complex parts: u = ui +iu2 
where Ui is a M- valued elements of Hs, etc. 
From (5.20), one computes 

(^l+i{u,0w-\{u,0lv + ---^ dP{0 = l-l{n,u)£ + --- . (5.21) 

Now it follows that E(u^) = E((u, C)vv) = ll^lll fo^' every u € 5£, by comparing the 
terms of (5.21) which are quadratic in u. Therefore, W : Us — > 5^ is an isometry, 
and (5.21) gives 

Ei\i,-iy\'')=E{{v,~Vy,0^) = \K~Vy\\l, (5.22) 

whence (5.18) follows from (5.8). Note that by comparing the linear terms, (5.21) 
implies E(l) = 1, so that P is a probability measure, and E((m,^)) = and 
E((w,^)^) = so that P is actually Gaussian. 

Finally, use polarization to compute 

{u,v)£ = ^ {\\u + v\\j - \\u-v\\l) 
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= - (E{\u + i\') -E(lfi-^I')) by (5.22) 

\u + i\^{^)-\u-v\^{^)dP{0 

This establishes (5.19) and, upon complexification, completes the proof. □ 

Remark 5.16. Observe that Theorem 5.15 was carried out for the free resistance, 
but all the arguments go through equally well for the wired resistance; note that 
is similarly negative semidefinite by Theorem 5.5 and [JP09c, Cor. 5.5]. Thus, 
there is a corresponding Wiener transform W : Ttn — )■ {S'g , P) defined by 

W:v^f, / - P^„,v and /(^ = (/, Ow- (5.23) 

Again, {fx}xeG° is a system of Gaussian random variables which gives the wired 
resistance distance by R^{x,y) — E((/a; — fy)^)- 

Remark 5.17. For u G Harm and ^ G S'^, let us abuse notation and write u for u 
so as to avoid unnecessary tildes. That is, u{^) := u{^) — (u,^)w. 

Remark 5.18. The polynomials are dense in L'^{S'g,¥): let (p{ti,t2, . . . ,tk) denote 
an ordinary polynomial in k variables. Then 

^{0 ■.^^{u,{0,U2{0,---Un{0) (5.24) 

is a polynomial on 5^ and 

Voly„:^{ip(uiiO,U2iO,---MO) ,deg{ip)<n, ■ e , ^ e S',} (5.25) 

is the collection of polynomials of degree at most n, and {Poly.^J'^^Q is an increas- 
ing family whose union is all of S'^. One can see that the monomials (u,Ow a-re 
in L'^{Sg,¥) as follows: compare like powers of u from either side of (5.21) to see 
that E((u, 0^+1) =Oand 

E((«,Ow)=/ l(^,OwrdP(0 = |?Tll"lir, (5.26) 

and then apply the Schwarz inequality. 

To see why the polynomials {7^oZy„}J^g should be dense in l?(S'^^V) observe 
that the sequence \Pvoiy^'^=Q of orthogonal projections increases to the identity, 
and therefore, {Ppoiy u} forms a martingale, for any u G (i.e., for any u G 

Denote the "multiple Wiener integral of degree n" by 

iJ„ := (dspan{(w, ] u G Hs}) G {{u,-)^ ■. k <n,u€ He}, 
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for each n > 1, and Ho CI for a vector 1 with ||1||2 = 1- Then we have an 
orthogonal decomposition of the Hilbert space 

oo 

L2(5^,P) = 0i?„. (5.27) 

Tl = 

See [HidSO, Thm. 4.1] for a more extensive discussion. A physicist would call 
(5.27) the Fock space representation of L'^{S'g,¥) with "vacuum vector" 1. Note 
that Hn has a natural (symmetric) tensor product structure: = Tif", the n- 
fold symmetric tensor product of Tig with itself. Observe that 1 is orthogonal to 
Hn and Harm, but is not the zero element of L'^{S'g,¥). 

Familiarity with these ideas is not necessary for the sequel, but the decom- 
position (5.27) is helpful for understanding two key things: 

(i) The Wiener isometry W : He — )■ L"^ {S'g , P) identifies He with the subspace 
Hi of L^(5£,P), in particular, L^(5£,P) is not isomorphic to "Hf. In fact, it 
is the second quantization of "He . 

(ii) The constant function 1 is an element of (5^ , P) but does not correspond 
to any element of He - In particular, 1 is not equivalent to in L'^{S£,P) (as 
it was in Hs). 

It is somewhat ironic that we began this story by removing the constants (via the 
introduction of He), only to reintroduce them with a certain amount of effort, 
much later. 

Recall that we began with a comparison of the Poisson boundary repre- 
sentation for bounded harmonic functions with the boundary sum representation 
recalled in Theorem 5.1: 

u{x)^l u{y)k{x,dy) O u{x) = ^ + u{o) . 



on 



bdG 



In this section, we replace the sum with an integral and complete the parallel. 

Corollary 5.19 (Boundary integral representation for harmonic functions). 

For any u € "Harm and with ~ PuarmVx, 



u{x)^ u{OK{OdV{0+u{o). (5.28) 
Proof. Starting with (2.6), compute 



u{x) - u{o) ^ {hx,u)£ ^ {u,hx)£ = / uhxdP, (5.29) 



where the last equality comes by substituting v — yci (5.19). It is shown in 
[JP09b, Lem. 2.24] that ^ h^. □ 

Remark 5.20 (A Hilbert space interpretation of bdG). In view of Corollary 5.19, 
we are now able to "catch" the boundary between and S'^: the boundary bdG 
may be thought of as (a presumably proper subset of) S'^/Hg. In parallel to 
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the construction of the Martin boundary, one expects that S'g/Hs is larger than 
necessary, and that P is probably supported on a much smaller set, comparable to 
the minimal Martin boundary; cf. [WocOO, Ch. IV]. Corollary 5.19 suggests that 
ii{x, c?^) := hx{£,)dP is the discrete analogue in Hg of the Poisson kernel k{x, dy), 
and comparison of (1.8) with (5.28) gives a way of understanding a boundary 
integral as a limit of Riemann sums: 

/ uKdP = \im u{x)^{x). (5.30) 

(We continue to omit the tildes as in Remark 5.17.) By a theorem of Nelson, P 
is fully supported on those functions which are Holder-continuous with exponent 
a = i, which we denote by Lip(i) C 5^; see [Nel64, Nel69]. Recall from [JP09c, 
Cor. 2.16] that Hg C Lip{^). Current research focuses on determining the precise 
relationship between bdG and these other spaces (Martin boundary, Lip{^)), and 
an explicit representation of bdG in terms of paths in G and/or cocycles. We 
expect that bd G will have applications in the analysis of self-similar fractals, by 
understanding the fractal as a boundary of a resistance network. 



6. Examples 

In this section, we introduce the most basic family of examples that illustrate our 
technical results and exhibit the properties (and support the types of functions) 
that we have discussed above. 

Example 6.1 (Geometric integer model). For a fixed constant c > 1, let (Z,c") 
denote the network with integers for vertices, and with geometrically increasing 
conductances defined by c„_i,„ = c™'*''^!"!'!""^!^ so that the network under con- 
sideration is 

3 2 2 3 4 

■ • ■ -2 -1 1 ^— 2 ^— 3 — 



Fix = 0. On this network, the energy kernel is given by 

'O, fc < 0, 

^^n(fc) = <!^^, l<fc<n,n>0, 



k > n, 



l-r"+i 



1-r 

and similarly for n < 0. Furthermore, the function 

u;„(n)=arH, a (6.1) 

defines a monopole, and h(n) = sgn(n)(l — Wo{n)) defines an element of T-Larm. 

Example 6.2 (Geometric half-integer model). It is also interesting to consider 
(Z+,c"), as this network supports a monopole, but has T-Larm — 0. 
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The monopole can be obtained by rescaling (6.1); just take a :— (iz^y There 
cannot be any nontrivial harmonic functions on this network by [JP09b, Lem. 5.5], 
which states that if ft, S 'Harm\{0} , then h has at least two different hmiting values 
at oo. That is, there exist infinite paths 71 = (xi, X2, ■ ■ ■) and 72 = (yi, y2j • • • ) 
with liiUj^oo K^j) 7^ limj^oo KVj)- 

For fc = 2, 3, . . . , the network (Z+, fc") can be thought of as the "projection" 
of the homogeneous tree of degree k {Tk, under a map which sends x to n S Z 
iff there are n edges between x and o. 

Remark 6.3. One can consider more general integer networks, and in this case. 
Harm ^ for (Z,c) iff "^c^y < 00. In this case. Harm is spanned by a single 
bounded function; details appear in [JP09b]. Networks of this form have been 
discussed elsewhere in the literature (for example, [KY89, Ex. 3.12, Ex. 4.9] and 
[KY84, Ex. 3.1, Ex. 3.2]), but the authors appear to assume that A is self-adjoint. 
This is not generally the case when c is unbounded; in fact, the Laplacian is not 
self-adjoint for Example 6.2 or Example 6.1; see [JP09d, §4.2] or [JP08, §13.4] for 
further discussion and the explicit computation of defect vectors. 

Example 6.4 (Star networks). Let (5„i, c") be a network constructed by conjoining 
m copies of (Z^., c") by identifying the origins of each; let o be the common origin. 

Recall from Theorem 3.10 that the boundary term is nontrivial precisely 
when bdG 7^ 0; the presence of a monopole indicates that bdG contains at least 
one point. If Harm ^ 0, then there are at least two boundary points; see [JP09b, 
Lem. 5.5] and Corollary 3.8. 

Example 6.4 shows how to construct a network which has a boundary with 
cardinality m. Note that these boundary points can be distinguished by monopoles, 
by constructing a monopole which is constant everywhere except on one branch. 

Example 6.5 (Networks of integer lattices). For d > 3, let {Z^^-,}^"^]^ be a collection 
of m copies of the rf-dimcnsional integer lattice Z'^ with edges between nearest 
neighbours, and let Ok denote the origin of Z^^,^. Let Z^ be the Cayley graph 
of the cyclic group of order m, and denote its elements by {l,2,...,m}. Now 
define (Z'^ ® Z^,!) by identifying Ok G Z^^.^ with k e Z„i, thus conjoining all 
the copies of Z''. Since Z"^ is transient, each copy Z^^,^ supports a monopole, and 
hence Harm has dimension m — 1 for this network. This is essentially a variation 
of Example 6.4 where (Z+, c") is replaced by Z*^. Note that this is not the same as 
the Cayley graph of the wreath product Z''' I Z„j, which is instead a Diestel-Leader 
graph; cf. [Woe05]. 

Example 6.6 (One-sided infinite ladder network). Consider two copies of the 
nearest-neighbour graph on the nonnegative integers Z"*" , one with vertices labelled 
by {xn}, and the other with vertices labelled by {j/n}. Fix two positive numbers 
a > 1 > /? > 0. In addition to the edges Cx^^x^^i — and Cy^^y^_^ = a", we also 
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R"\x,o) = &{f,) = \ 



R"(x,o) = S{h,) = : 



Figure 3. The reproducing kernel on the tree with c = 1. For a vertex 
X which is adjacent to the origin o, this figure illustrates the elements Vx, 
fx = Pj-inVx, and hx = P-HarmUx; See Example 6.7. 



add "rungs" to the ladder by defining Cx 



X2 



X3 



13 






a , 







Vi 



y2 



2/3 



Vn 



(6.2) 



This network was suggested to us by Agelos Georgakopoulos. In [JP09a], we show 
that this example is a one-ended network with nontrivial "Harm, by explicitly 
constructing a formula for a harmonic function of finite energy on this network. 

Example 6.7 (The reproducing kernel on the tree). Let (T, 1) be the binary tree 
network as in the top of Figure 2 with constant conductance c = 1. Figure 3 
depicts the embedded image of a vertex v^, as well as its decomposition in terms 
of J-ln and T-Larm. We have chosen x to be adjacent to the origin o; the binary 
label of this vertex would be xi. 

In Figure 3, numbers indicate the value of the function at that vertex; artistic 
liberties have been taken. If vertices s and t are the same distance from o, then 
|/x(s)| = I/kI^)! similarly for h^. Note that hx provides an example of a 
nonconstant harmonic function in %£. It is easy to see that Ymiz^±ao hx{z) = |±i, 
whence h^- is bounded. 
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We can use of Figure 3 to describe an infinite forest of mutually orthogonal 
harmonic functions on the binary tree. Let z G 7~ be represented by a finite binary 
sequence: the root o corresponds to the empty sequence 0, and the two vertices 
connected to it are and 1. The neighbours of are 0, 00 and 01; the neighbours 
of Of are 0, OfO, and Off, etc. Define a mapping ip^ : T T by prepending, 
i.e., ipz{x) = zx. This has the effect of "rigidly" translating the the tree so that 
the image lies on the subtree with root z. Then hz '■= h^oipz is harmonic and is 
supported only on the subtree with root z. The supports of h^^ and h^^ intersect if 
and only if Im(c/J2. ) C Im((^2.). For concreteness, suppose it is 1tii{lPz-^) Q lu\{ipz^). 
ff they are equal, it is because zi — zi and we don't care. Otherwise, compute the 
dissipation of the induced currents 



Note that dh^^ix^y) always has the same sign on the subtree with root z\ ^ o, 
but dhz-^{x,y) appears in the dissipation sum positively signed with the same 
multiplicity as it appears negatively signed. Consequently, all terms cancel and 
= {dhzi,dhz2)D = {hzi,hz2)£ shows hz^ -L hz^- 

This family of harmonic functions can be heuristically described by analogy 
with flaar wavelets.^ Consider the boundary of the tree as a copy of the unit 
interval with as the basic Haar mother wavelet; via the "shadow" cast by 
lim„^±oo hx{xn) = ±f (this can be formalized in terms of cocycles). Then hz is 
a Haar wavelet localized to the subinterval of the support of its shadow, etc. Of 
course, this heuristic is a bit misleading, since the boundary is actually isomorphic 
to {0, 1}^ with its natural cylinder-set topology. 
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